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Using the adaptive time-dependent density matrix renormalization group method, we numerically 
investigate the expansion dynamics of bosons in a one-dimensional hard-core boson model with 
three-body interactions. It is found that the bosons expand ballistically with weak interaction, 
which are obtained by local density and the radius R„. It is shown that the expansion velocity 
V, obtained from R„ — Vt, is dependent on the number of bosons. As a prominent result, the 
expansion velocity decreases with the enhancement of three-body interaction. We further study the 
dynamics of the system, which quenches from the ground state with two-thirds filling, the results 
indicate the expansion is also ballistic in the gapless phase regime. It could help us detect the phase 
transition in the system. 
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INTRODUCTION 

Recently, the understanding of nonequilibrium dynam¬ 
ics of strongly correlated many-body systems poses one of 
the most challenging problems for both theoretical and 
experimental physics Researches in the nonequilib¬ 
rium properties of strongly correlated many-body sys¬ 
tems have emerged into a dynamic and active held, 
driven by the possibility to address questions such as 
thermalization EHl] and particle transportation [SHia in 
clean, well-controlled, and isolated systems. In the par¬ 
ticle transport cases, there are two prototypical trans¬ 
port mechanisms in classical physics: ballistic and dif¬ 
fusive transports. The ballistic and diffusive transports 
are characterized by non-decaying currents and decay¬ 
ing currents respectively. In the ballistic systems there 
is absence of friction. However, in the diffusive systems, 
there have frequent diffractive collisions, which drive a 
local thermalization. Many studies focus on qualitative 
questions such as whether transport is ballistic or rather 
diffusive in microscopic models of strongly interacting 
systems 1 141119] . In experiment, the measurements of lo¬ 
cal occupancy dynamics can be realized even the densities 
of initial states are than one[5D]. The dynamics exhibit 
clear signatures of quantum distillation and conhnement 
of vacancies in the doublon sea, which are significantly 
interesting. Other typical examples investigated numer¬ 
ically are the expansion of initially localized ultra-cold 
bosons in homogeneous one- and two-dimensional opti¬ 
cal lattices. It is found that both dimensionality and 
interaction strength crucially influence these nonequi¬ 
librium dynamics, which have also been confirmed in 
experiments |21| . 

Previous results for nonequilibrium dynamics strongly 
correlated many-body systems are almost with the dom¬ 
inant two-body interactions, because relatively small 
multi-body interactions can only provide tiny corrections. 


Recently, the cold atom in optical lattice gives us a great 
platform to realize three-body interactions [52H21]- It 
is shown that the three-body interactions can be dom¬ 
inated, and the two-body interactions can be indepen¬ 
dently controlled and even switched off by driving mi¬ 
crowave fields|22]. The system with multi-body inter¬ 
actions can induce many exotic phenomena. It would 
be interesting to investigate nonequilibrium dynamics in 
the system with multi-body interactions, such as sudden 
expansion of Mott insulators (MI) in a one-dimensional 
hard-core boson model with three-body interactions. 


The paper is organized in the following ways. We de¬ 
fine the model Hamiltonian, and the observable is also 
provided. The evidence of the difference between bal¬ 
listic and diffusive expansion dynamics, which can be 
obtained by the behaviours of the observable, is dis¬ 
played. Then we focus on the sudden expansion in a one¬ 
dimensional hard-core boson model with the three-body 
interaction. The results presented in the main text are 
obtained for the expansion from the box trap. The hard¬ 
core bosons that expand from the product of local MI 
states is probed. Furthermore, we investigate the hard¬ 
core bosons expansion from the entanglement states. A 
discussion follows in the last section. 


As a result, when the hard-core bosons expand from 
the product MI states, we obtain the bosons expansion 
ballistically with weak interaction. The expansion veloc¬ 
ity decreases with the increase of three-body interaction. 
Moreover, we study the dynamics in the system from the 
ground state with two-thirds filling, and our results in¬ 
dicate that the expansion is also ballistic in the gapless 
regime. 
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FIG. 1. (Color online) Typical contour plot of the local density {ni{t)) as function of position and time during the expansion 
from a MI(N=20,L=180) with different three-body interaction (a)VF = 0,(b)IF = 1 and {c)W = 3. The slanted lines in (a) and 
(b) indicate the speed 2J. It is noted that the strength of (a),(b) and (c) has the same colour-bar. 


MODEL AND MEASUREMENTS 


The Hamiltonian we consider is one-dimensional hard¬ 
core boson model with three-body interactions, and is 
given by 

H = -jy^{blbi+i + bibl_^_^) + Wy^nj-iniTii+i, (1) 

i i 

where is creation (annihilation) operator of hard¬ 

core boson at site i: (6j)^ = (6^)^ = 0, and rii = b\bi 
is local density operator. The parameter J is the hop¬ 
ping interaction and chosen as the unit of energy in the 
paper, and only the leading three-body interactions with 
strength W is considered. Open boundary conditions are 
imposed in the system. 

As we known, there are two prototypical transport 
mechanisms in classical physics: ballistic and diffusive 
transports. The two kinds of transports can be distin¬ 
guished by using the time dependent radius of the density 
distribution, which is defined as 








Rl(t = 0), 


( 2 ) 


where the parameter N indicates the number of the par¬ 
ticles and L indicates the length of lattice sites. The pa¬ 
rameter iQ represents the center of mass ig = Ljl+Xj^. A 
ballistic expansion will lead to the well-known i?„(t) oc t 
behavior, while a diffusive expansion with a fixed diffu¬ 
sion constant, will lead to the well-known Rn{t) oc '/i 
behavior. The radius has the ability to detect trans¬ 
port whether if is ballistic or diffusive. These have also 
been verified for spin and energy dynamics in the spin- 
1/2 XXZ chain 0115]. 

Several powerful methods have been employed to study 
the expansion dynamics in one dimension, including the 
time-dependent density matrix renormalization group 



FIG. 2. (Golor online) The radius Rn{t) is plotted as function 
of time for different W with N = 20. 


method (t-DMR,G) |25l[2i^ . The method is also applied in 
the paper. In our simulations, the Trotter slicing Jt = 0.1 
is used, and the t-DMRG codes with double precision are 
performed with a truncated Hilbert space of m = 400. 
The time that can be simulated in the system is deter¬ 
mined by the entanglement entropy. 


EXPANSION FROM MI STATE 


For simplicity, we hrst discuss the idealized case that 
all the particles are in a box-like trap. It is not difficult 
to achieve in experiments, just like by using the confin¬ 
ing potential m- It is important to emphasize that the 
confining potential would influence the three-body inter¬ 
actions. After relaxation of the confining potential, its 
effects on the three-body interaction will disappear. As 
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FIG. 3. (Color online)The total particle current Ji /2 vs time 
for different three-body interaction W. 



FIG. 4. (Color online) The expansion velocity V vs W with 
different N. The inset: Finite-size scaling of V{N) with W = 
l(o) and W — 2(n). The lines are the fit lines. 

is known, the theoretical investigation for two-body in¬ 
teractions in such a system is very interesting, and has 
been investigated in Refs. [33 [31]. Here, we will focus our 
attention on the effects of three-body interactions. The 
initial states are that all the particles are in a box trap 
states {Lbox is the length of the box) and the box trap is 
full of with particles, and the initial state is a product of 
local MI state in the box trap, and is given by 

\'^^n^Ual)= n ( 3 ) 

The parameter N indicates the number of the particles 
N = Lbox = *2 — H + 1, and this state has been realized 
in a recent experiment |2Ij. 

A typical example for the time evolution of the density 
{ni{t)) is shown in Fig. [^for the expansion from a MI 
state with different W. In Fig. 0a) , it is found that the 


MI melts on a time scale of tmeit — Aboa;/2/(2J) = 5 for 
W = 0. For the case of IT = 0, the Hamiltonian ([^ can 
be diagonalized by using Jordan-Wigner transformation 
and Fourier transformation, which is given by 

H = '^€kclck, (4) 

k 

where 

Cfc = —2Jcos{k), k = Itt/{L + l),l = 1,2,3 ■■■ L. (5) 

Therefore 2J is the largest possible velocity. In Fig. [T] 
(a) and (b) two lines are parallel to the outer rays visible 
in the figure, i.e., the fastest propagating particles, and 
it indicates an excitation spreading out from the center 
box with the same group velocity 2 J. When IF = I, tmeit 
increases, as shown in Fig. 0b)- For t > tmeit, two par¬ 
ticle clouds form that propagate into opposite directions, 
visible as two intense jets, as is expected for ballistic dy¬ 
namics in one dimension system [TU [T^. When IF = 3, 
we cannot find tmeit in Fig. [^(c). With strong three- 
body interaction, the particles are almost all confined in 
the box trap. There is no possibility of compensating 
the loss in interaction energy with kinetic energy in the 
system due to energy conservation. 

In order to clearly distinguish whether the expansion 
is ballistic or not, we also investigate the radius behav¬ 
iors. In Fig. we show the radius R„(t) for differ¬ 
ent magnitude of three-body interactions. Clearly when 
IF = 0,1.0,1.6, the radius i?„(f) ^ t within short and 
intermediate time. In the time region, the local densities 
are large, and the expansion could be considered ballis¬ 
tic. For IF = 2.0, 3.0, in the initial time it does not follow 
Rn{t) ^ t, and after that the radius follows Rn{t) ~ t. It 
seems that the expansion is ballistic obtained by the ra¬ 
dius. But in general, ballistic dynamics in strongly inter¬ 
acting dynamics is connected to integrability. The hard¬ 
core bosons without any additional three-body interac¬ 
tions expand ballistically due to their integrability and 
the mapping to non-interacting fermions that preserves 
density. The three-body interactions break integrability 
and the dynamics would be expected as diffusion. With 
strong three-body interactions, that the particles are al¬ 
most all confined in the box trap also suggests that the 
transport would be diffusive [33] . However, the sudden 
expansion of initially trapped particles into an empty lat¬ 
tice is more complicated, this can not provide sufficient 
evidence for diffusion. To be sure, the intermediate time 
regime in which i?„ is not proportional with time could 
be an indicator of non-ballistic dynamics. 

In order to check the results for IF = 2.0,3.0, we inves¬ 
tigate the time dependence of the total particle current 
in each half of the system, which is dehned as 

JL /2 = -i X! J{4bi+i- h.c.). (6) 

i>L/2 

It is noted that ballistic dynamics are probed into the 
systems with open boundary conditions due to the exis¬ 
tence of globally conserved currents[33|- The results of 
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Jl/2 for W = 2,2, are shown in Fig. The parameter 
Jl/2 takes a constant value for W =^,3, which is the 
benchmark feature of ballistic dynamics [T5]. However, 
the behaviour of Jl/2 does not imply that the dynamics 
will be ballistic. There is a subtle solution of the equation 
that the radius increases linearly with time (if properly 
defined in Eq. and the total particle current in each 
half of the system takes a constant. The phenomena gen¬ 
erally occur in the dilute limit. With strong interaction 
W = 2,3, the particles are almost all confined in the box 
trap, and only a few particles transfer. This case could be 
considered the dilute case, and the dynamics will expand 
ballistically at a very long time, irrespective of whether 
it is ballistic or diffusive mass transport. 

We extract the expansion velocity V by fitting the 
tDMRG data , i.e., the slope of curves such as the ones 
shown in Fig. with Curve Fitting toolbox in Matlab. 
The corresponding radial velocity V is defined through 
the reduced radius V = dRn{t)/dt. As is known, there 
are two main sources of errors in the adaptive t-DMRG: 
the Trotter error due to the Trotter decomposition and 
the DMRG truncation error due to the representation 
of the time-evolving quantum state in reduced Hilbert 
spaces. For small times, the Trotter error plays a major 
role. For long times, the DMRG truncation error will 
be dominated[29l [34]- We obtain V with the date 
in time range 2 < t < 10, where both Trotter error and 
truncation error are small. Results for selected values 
of W are collected in the main panel of Fig. For 
the expansion from the MI we obtain V = V2J for any 
A' > 0 and W = Q [IH]. The expansion velocity also 
decreases with the increase of W for the same A. This 
results in an increase of interaction energy and therefore 
a decrease in kinetic energy. The velocity V, which em¬ 
bodies the behaviours of kinetic energy, decreases when 
the interaction increases. The expansion velocity is much 
smaller than the largest possible velocity 2 J and is always 
very different from characteristic velocities of the initial 
state. Fig. also indicates that the expansion velocity 
decreases while A increases, and the finite-size effect is 
shown in the inset. For above case, the length of system 
does not affect the propagation noticeably until the par¬ 
ticles wave front reaches the boundary. We can fit the 
locations of maximums by the formula 

V(N)^V, + aN-\ (7) 

where a is a size-independent constant and A the number 
of bosons[T7]- We obtain that Wi = 0.95 for IF = 1 and 
Vc 2 = 0.30 for IF = 2, see Fig. inset. The results 
A —7^ oo for different IF is also shown in Fig The 
above results about the expansion velocity can provide 
some information about the initial state. 


EXPANSION FROM THE GRAND-STATE 

In this section, the initial state can be realized by find¬ 
ing the ground state of a system with (A — 2) and then 



FIG. 5. (Color online) The radius Rn{t) is plotted as function 
of time for different W with (n) = N/L — 2/3. 



FIG. 6. (Color online) The expansion velocity F is plotted 
as function of the three-body interaction IF with L = 180. 
Inset: Finite-size scaling of critical Fc(A). The lines are the 
fit lines. 


applying two operator flips b^L/ 2 ^^L/ 2 +i- time evo¬ 

lution is then performed at the two-thirds filling (n) = 
A/L = 2/3. It is noted that with this filling, the quantum 
phase transition from the superfluid to solid phases in the 
system is located at Wc — 2.80±0.15, which are obtained 
by the structure factors and bond order Darameter|35) . 
In order to find whether the expansion is ballistic or dif¬ 
fusive, we also investigate the radius behavior. In Fig. 

we display the radius Rn(t) for various three-body in¬ 
teractions for (n) = 2/3. Glearly when W = 0,1.0,2.0, 
the radius is proportional with timei?„(t) ^ t. It seems 
that the expansion is ballistic. It is noted that the radius 
obtained by Eq. ([^ does not follow oc t or (x Vl when 
W > 2.8 due to the strong oscillations of {ni{t)) [14]. We 
also adopt the method of Eq. (13) in Ref. jHj, the os¬ 
cillations also exist. We cannot conclude the transport is 
ballistic or diffusive by the radius. We extract the expan- 
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sion velocity V by fitting the tDMRG data. We obtain 
V with the data in time scale 2 < t < 10, which 

the Trotter error and truncation error are both small. 
The velocity for different values of W are collected in the 
main panel of Fig. It is seen that the velocity has 
a peak which is located at W ~ 2.5. The initial states 
would be influenced by the finite size effect. The loca¬ 
tion of peak may change for different system size. So 
we investigate the location of peak with different system 
sizes. It is found that the locations of expansion velocity 
maximums follow the formula 

W(L) ^ Wc + cL-\ (8) 

where c is a size-independent constant and L the sys¬ 
tem size, seen in Fig. [^inset. It is found that the critical 
point Wc — 2.74±0.1, which is close to the the superfluid 
to solid phases quantum phase transition point. The dy¬ 
namics would be diffusive when W > |36If55] . where 

the system is in the gapped phase [3H1. It means that 
the peak of velocity will move to the boundary between 
two kinds of transport behaviors, as the system size ap¬ 
proaches to infinite. 

DISCUSSION 

In the paper, we study the expansion Mott insulator 
in a one-dimensional hard-core boson model with three- 
body interactions by using the adaptive time-dependent 
density matrix renormalization group method. We ob¬ 
tain that the bosons expand ballistically with weak inter¬ 
action by studying the dynamics of local density and the 
radius i?„. It is also found that the expansion velocity V, 
is dependent on the number of bosons in the system, and 


its measurement can provide information about the ini¬ 
tial state. As a result, the expansion velocity decreases 
when the strength of three-body interaction increases. 
Moreover, we study the dynamics of the system from 
the state with the two-thirds filling. Our results indicate 
that the expansion also is ballistic in the gapless phase. 
In the case of the two-thirds filling, the superfluid to solid 
phases quantum phase transition point can be identified 
by the peak of the expansion velocity. It would be in¬ 
teresting to test our predictions in experiments. With 
strong interaction, the dynamics could not be ballistic. 
It would be significance work to get the expansion dy¬ 
namics by other parameters. 
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